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Abstract 

' With the latest LHC available results, we consider the generic constraints on massive graviton. 

Both dijet and dilepton resonance searches are used. The limits on parameter space can be applied 

(N 

to many models. As an illustration, we show the constraints for Randall-Sundrum (RS) model. 
Implications on massive graviton and the coupling strength are discussed. For k/M p i = 0.1, 
Mq < 2.2 TeV region is excluded at 95% confidence level. We also present some interesting 
^ <— | implications on the RS radion with respect to the 125 GeV excess at the LHC. For k/M p i = 0.1, 

A<£ < 13.8 TeV is excluded where A^ is the scale to charactarize the interaction strengh of radion. 
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I. INTRODUCTION 



The standard model (SM) with gauge group SU(3) xSU(2) xU(l) has been very successful 
in explanation for experimental results. Although the exact mechanism for electroweak 
symmetry breaking is still unclear, the recent hint of 125 GeV excess jl, 2] may suggest the 
higgs mechanism. Then the hierarchy problem in SM, electroweak breaking scale being so 
smaller than Planck scale, still exists and motivates new physical ideas. Among these ideas 
to solve the hierachy problem, extra-dimension {3, 4] is one of the simplest ways . 

Numerous models with extra-dimensions describing possible new physics beyond SM 
predict massive gravitons {|3, 6]. The masses of these particles are usually at the TeV scale 
if these models were trying to solve the hierarchy problem. Then these massive gravitons 
can be produced at the large hadron collider (LHC). In some cases, the cross section is large 
enough so that exclusion limit or discovery can be reached for the considerred models. 

Extra-dimensional models have many phenomenolgical consequences UQ. Before LHC 
era, flavour physics and electroweak precision observables have already give some limits. In 
our study, we shall discuss the constraints from direct searches at the LHC. Other than 
specifying on a single model, we will study the general properties of massive graviton. The 
discussions and constraints are applicable to a wide class of models. As an illustration, we 
show the constraints for Randall- Sundrum (RS) model. 

The paper is organized as follows. In section |TT1 we introduce the general framework 
for later discussion and the necessary ingredients for searches or exclusions for massive 
graviton at the LHC. In section IHIl we use the latest CMS and ATLAS data to constrain 
the parameters for the lightest massive graviton. Both dijet and dilepton final states are 
considered. In section IIVt we show that constraints has interesting implications in a specific 
and popular model, warped extra dimension. Finally, summary and conclusion are given. 



II. GENERIC FRAMEWORK ON MASSIVE GRAVITON 

Massive gravtions exist in various models usually with extra dimensions. The compacifi- 
cation of extra-dimensional leads to Kaluza-Klein towers of particle spectrums. The lowest 
states of these towers usually are the SM particles, and the higher states represent these new 
heavy particles. With exact mass depending on the details of models, massive Kaluza-Klein 
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particles are roughly at TeV scale if heirarachy problem is solved. Among them, massive 
gravitons serve as indispencible ingredients. 



A. Interactions with SM 



The field that describes a spin-2 particle is a tensor, h^. In a general effective theory 
with h^v, all gauge invariant terms, both renormalizable and non-renormalizable, should 
be written down in the lagrangian. However, such a theory will have to too many free 
parameters to be considered. Thus, in this paper, we should assume a minimal setup that 
massive graviton shall couple to standard model particles with universal form as same as 
the massless one, except with different strength and a non-zero mass. This universality is 
quite general as long as massive graviton is derived from the space-time metric. 

For direct production at the LHC, we only need the lowest interaction term, 

Ant = ~^T a "h% - ^f^hap, (2.1) 

here and after, we will use hP and h^ v for massless and massive graviton, respectively. 
T a/3 is the energy-momemtum tensor of SM field. And the interactions between massive 
graviton and SM particle are solely determined by Ac- In general, there could be more than 
one massive graviton. For simplicity in this work, we will only concentrate on and refer to 
the lowest massive graviton, h^^jS not stated explicitly. Early discussions on searches for 
massive graviotn are referred to 13K17]. 



B. Production at the LHC 



Based on the interaction term, it is immediately realized that the topology for production 
of massive graviton is Drell-Yan like process. At the LHC with proton-proton collision, the 
two main channels to produce massive graviton are gluon-gloun fusion and quark-antiquark 
annihilation processes, Fig. (JT|). Vector boson fusion processes or other assiociated produc- 
tion are neglected for leading-order approximation throughout this work. 

For hadronic collison, the total cross section is the convolution of the parton distribution 
fucntions (PDFs) with the partonic cross section, 

a = J dxidx 2 f qi (x 1 ,fi F ) f q2 (x 2 , fJ, F ) a (<?i<?2 -»■ G*; s) , (2.2) 
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FIG. 1. Feynman diagram for the massive graviton production at the LHC. The produced gravitons 
then decay to SM particles. 




FIG. 2. Colorred map of cross section (pb) for massive graviton production at the LHC with 
y/s = 7 TeV (left) and y/s = 8 TeV (right), respectively. 



where f q (x,/j,F) is the PDF for a parton g(quark and gluon) with momentum fraction x at 
the factorization scale fip, a is the partonic cross section with the initial two partons of 
momentum fraction x% and X2, respectively, and s = X\X2S. 



We use Madgraph 5 (l8|, ll9| with fi F = M G and CTEQ6L1 20] PDF set, and show in 
Fig. (J2]) colorred maps of cross section with function as mass Mq and the scale at the 
LHC with both y/s = 7 TeV (Left figure) and y/s = 8 TeV. The cross section smaller than 
1 fb is shown in blue region at the right top corners. These two figures roughly show the 
relations between cross section and parameters. Compared to the y/s = 7 TeV case, the 
colored map for y/s = 8 TeV is shifted towards larger Ac and Mq. 



C. Decay width and branching ratio 



For direct search or exclusion of the massive graviton, we need to know the branching 
ratio of its decay channels. If massive graviton lies in the TeV mass regions, its decay 
products, standard model particles, are then on-shell. So for leading order consideration 
we only include the two-body decay channels. The individual decay rate to two final SM 
particles is listed in the appendix. In Fig. [3j we show the branching ratios for the main 
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FIG. 3. Branching ratios, where q = u,d, c, s, b and I = e, \i. We have Br(G — > 11) ~ Br(G — > TT) 
as shown above. 



decay channels. Here, we neglected the three-body decay with a off-shell W*, Z* or t at the 
low mass range. As shown in the figure, as the mass of graviton goes large compared with 
2mt, branching ratios are almost fixed. Now it is easy to get the total decay width for large 
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In Table. [U we show branching ratios of the main decay channels for later use. The events 
at partonic level at the LHC then is proportional ot cr(pp — > G) x Br(G — > //) for narrow 
width approximation. 



III. CONSTRAINTS FROM LHC DIRECT SEARCHES 



The massive graviton couples to standard model particles through the energy-momentum 
tensor in the linear theory and the couplings are suppressed by a factor Aq- If the mass is at 
the TeV scale and the factor is not too large, then enough massive graviton can be produced 
at the LHC. Searches for the decay products of the graviton can be used to discover or 
constrain the parameter space of the model. Based on the final states, there are various 
search channels and if no excess is observed, each channel can give a constraint. In the 
section, we focus on the dijet and dilepton final states only. 



A. Dijet Constraint 

Dijets consist of quark jets and gluon jets. The shape of a gluon jet is wider than that 
of a quark one because gluon's effective coupling Caois{Ca = 3) is larger than quark's 
CfC(s{Cf — 4/3), then gluon jets are more likely to radiate. Modern detectors have the 



power to distinguish them. CMS [21] has presented limits for three kinds of dijet resonances, 
gluon-gluon, quark-quark( quark for q or q) and gluon-quark. Because of the larger back- 
ground from gluon jets, quark-quark dijet limit is the most stringent one among three . 
Then, we shall only consider the constraint from quark-quark dijet spectrum. 

From the branching ratio in Table. [H we know that about 1/3 of the produced massive 
graviton will decay to quark-antiquark pairs (qq, q = u,d, s,c,b). If the production rate is 
large, event distribution on dijet invariant mass rrijj shall show an additonal peak at the 
graviton mass in the smooth QCD dijet background. If no excess is observed, constraints 
can be put on production rate using statistics(see appendix). So far, the latest published 



model-independent result on resonance search by dijet search is from CMS with 1 fb 1 21 ]. 



Since no excess is observed yet, we sha 



1 use this result to constrain our parameter space. 



The analysis of event samples in |21( relies on the following selection rules, 

p T > 10 GeV, \rj\ < 2.5, |A??| < 2.5 and m# > 838 GeV, (3.1) 



TABLE I. Branching ratio for the massive graviton to SM particle with I = e, \i and q = u,d, s, c, b. 



channels 


HH 


99 


77 


W + W~ 


ZZ 


tt 


qq 


i+r 


Br 


2/293 


96/293 


12/293 


24/293 


12/293 


16/293 


90/293 


12/293 



6 



where the definitions are 



= a/ (E 1 + E 2 ) 2 - (pi + p 2 ) 2 , Pt= \ pl+pl=P sin 9 



6 A 

T) = - In tan -, Arj — r} x - rj 2 . 

The —nation S ca,e * L to „ = PrJP d CTEQ6L1 p art o„ dlstrlbutl on functl on S U 
are used. A K-factor of 1.33 was used in 



211]. 



For mass graviton, the signal acceptance A is about 0.72 for the events selected that 
satisfying the above kinematics requirements. This factor is nearly constant for large mass 
of G. Without seeing any excess of dijet mass spectrum, an upper limit at the 95% confidence 
level is obtained for a x Br x A, the products of cross section, branching ratio and events 
acceptance. Since both Br and A are known, the limit then is translated to constraint on 
the A G and M G . 

In Fig. HI we show the limit as a solid line on the cross section contour of and Mq. 
Regions at the left-handed side of the line is excluded at 95% confidence level. As is displayed 
in the figure, for Mq less than 1 TeV, A^ has to be larger than 2.5 TeV. As the mass get 
smaller, Aq needs to larger to accommadate. This is reasonable and intuitive otherwise we 
shall have seen a resonance around Mq- 



B. Dilepton Constraint 

The dilepton searches are involved with dielectron and dimuon final states where the 
main SM background is from Z/^* decay. For massive graviton, although the branching 
ratio to dilepton (/ = e,fi) is much smaller than that of dijet (7.5 times smaller exactly), 
the constraint or discovery potential is better due to the well-known and lower background. 



The event selection follows 22[ where for dielectron 

E T = sjp^ + m 2 > 25 GeV, \r)\ < 2.47 with excluding 1.37 < < |1.52|, (3.2) 

and for dimuon px > 25 GeV. The renormalization scale is set to [i = v§, the K-factor for 
dilepton final state varies between 1.6 to 1.8, depending on the graviton mass and A^. In 



practice, 1.75 is used for Mq > 750 GeV [22[. With these kinematic requirements, the total 
signal acceptance A is 72% for dielectron and 47% for dimuon. 

As shown in Fig. [4] as the dot-dashed line, the dilepton final state can give more stringent 
constraint due to both larger luminosity and better discriminating power, compared with 
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FIG. 4. No obsearvation of dijet and dilepton event excess can put constraints on the parameters. 
Cross section is in pb. The solid line indicates the constraints from CMS with 1 fb" 1 data [21] and 
the dot-dashed line shows the limit put by ATLAS with 5 fb -1 22]. 



dijet limit with only 1 fb data. The limit on cross section is almost constant 2 TeV and 3 TeV 
because the constraint mainly comes from the observed event H] in a single bin [1200,3000] 
GeV of [22]. 



IV. CONSTRAINT ON WARPED EXTRA DIMENSION 
A. Theory Overview 

One of the most popular extra-dimensional models is Randall- Sundrum model. The 

fl 

physics behind RS model lies in the following geometry for the warped space-time [4J, 

ds 2 = e - 2kT{x)M [r]^ + G^x^dx^dx" + T 2 {x)d<p 2 , (4.1) 

where T(x) is referred to as the modulus field, G^x) as graviton and k is a scale of the order 
of the (reduced) Planck scale M p i. To explain the hierarchy problem, the compactification 
radius or the vacuum expectation value(vev) of the modulus field, r c = (T(x)), is required 
to satisfy the relation kr c ~ 12. 
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The action that determines the above geometry is 

S = -Ml J d 5 xA^R {5) , (4.2) 
where the 5D reduced Planck scale M* is related with the 4D (reduced) Planck scale M p i, 

Ml = Ml / e~ 2k ^dy = ^ (l - e~ 2kr ^) . (4.3) 

J y=—r c n ™ 

In the original RS model, SM particles are confined to the brane and only graviton 
can propogate in the bulk. Then a massive Kaluza-Klein graviton tower exists besides the 
massless graviton, 

il4 n) = kx n e~ kr ^ = x n ^-A G , Aq = M pl e~ k ^\ J^O = 0, (4.4) 

where J\ is the Bessel function, and x\ ~ 3.8317, x 2 — 7.02, x 3 ~ 10.17, and x 4 ~ 13.32. 
We shall only consider the effect of the lowest state M G = Mq. Both massless and massive 
gravitons can couple to standard model particles, 



C 



M~T K ?~~A^ ^=iKp- (4.5) 



Mq = Xijj-Ac implys that the larger Mq is, the larger ratio k/M p i for fixed A G . 

pi 

Using the relation between the scale Ac and the mass Mq, we have the decay width for 
massive gravtion in RS model, 

293M r 293x?M G / k V „ / k \ 2 . . 

T G ~ f- = ~ 1.425M G . (4.6) 

G 960ttA g 960tt \M pl ) \M P J K J 

In the RS original paper, was assumed to be less than 1. Most discussions lie in 

0-01 < < 1 for theoretical and experimental studys. The estimation goes as follows. 

The 5D curvature scalar is R§ = —20k 2 , and requiring \R\ ~ Ml with Eq. 14.31 gives 

3 

20A; 2 ~ (kM 2 ^ =► ±- ~ (1 ) 4 ~ 0.1. 



However, it was arg ued in [23|, |2J, |33j that R 5 should be compared with A 2 (A = 24 1 / 3 ttM* 



is the energy scale at which the 5D gravity theory becomes strongly coupled), giving 

3 

20k 2 < A 2 =>- — — < | 24 2/3 vr 2 x — V ~ 2.88. 



M pl V 20 

When k/M p i is large, the decay width in Eq. 14.61 shows that it can be comparable with or 
even large than its mass. This can happen for strong interaction. For example, in PDG {25 1 
m/o{600) = (400 - 1200) MeV, but its width is in the range, T = (600 - 1000) MeV. And 



m p(770) = 775.49 MeV with width T = 149.1 MeV 
considered in the following section. 



The finite decay width effect will be 



B. Limits on RS graviton 
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FIG. 5. Constraint on the parameters Mq and Aq with dijet events at CMS 1 fb 21] 



G ->l + l",s =7 TeV 



ATLAS 5 fb" 
k/Mp,=0.1 
k/M pl =0.2 
k/Mp,=0.3 
k/Mpi=0.4 




1000 1500 2000 2500 3000 3500 
M G (GeV) 



FIG. 6. Constraint on the parameters with dilepton events Mq and Aq, for various k/M p i. The 



box points are extracted from 
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For the lightest massive graviton in RS model, we have the relation Eq. I4.4[ Mq = 
Xij^Aq. This feature shows that for fixed Mq, k/M p i can effectively discribe the interaction 
strength. Larger k/M p i means smaller and then stronger interaction. Both Tevatron and 
LHC has set some exclusion limit for 0.01 < -M- < 0.1. From the latest dilepton search 



result [22j, the k/M pl = 0.1 gives Mq > 2.16 TeV, implying A G > 5.64 TeV. 

In the following discussion, we shall extend the limit to k/M p i > 0.1 region. As shown in 
Fig. [5] and El both the dijet and dilepton exclusion limits depend on the k/M p i. In the dijet 
case, Mq < 1 TeV is excluded for k/M p i = 0.1 and the limit turns higher mass for larger 
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k/Mpi. In the dilepton figure, Mq < 2.2 TeV is excluded for k/M p i = 0.1, and again larger 
k/Mpi gives even more stringent constraints. 

For Drell-Yan like s-channel production followed by immediate decay, there is a propoga- 
tor of Breit-Wigner like 

1 

s- MI + iM g Tg 

where the width of the massive particle Tq has been included. Usually, the width is very 
small for weak interactions, the narrow width approximation can be used and the final cross 
section is product of cross sction for massive particle and branching ratio to final states. For 
large decay width, the above full Breit-Wigner propogator is needed. 

How the events are distrubted with respect to the van depends on PDfs f q (x,/1f), Mq 
and k/Mpi. In the Fig. [7J we show several cases with different Mq and k/M p i. When 
k/Mpi < 0.3, a clear resonance is still visible due to Tq < 0.2Mg for k/M p i < 0.38. When 
k/Mpi is large, the decay width goes large and the resonance get broadened. 

The limits are shown in Fig. [HJ For Mq > 3.5 TeV, the constraint is allmost fixed for 
large k/M p i. This is simply due to k/M p i cancellation between the coupling and decay width. 
And the processes can be discribed by a four-fermion contact interaction. 




1000 1500 2000 2500 3000 3500 4000 1000 1500 2000 2500 3000 3500 4000 

m„(GeV) m„(GeV) 



FIG. 7. Differential cross section as distributions with invarirant mass rrijj. Two cases with 
Mq = 2.5 TeV(left) and M G = 3 TeV(right) are shown. 



11 



M G (GeV) 

FIG. 8. Conservative limit on k/M p i and Mg, where the shadowed region is excluded at least at 
95% confidence level and the region between two solid lines indicates the effect of 10% uncertainty. 

C. Implications for RS Radion 

Randall-Sundrum scenario 4J is proposed to solve the hierarchy problem of the standard 
model. In the original version, no mechanism for dynamical origin of the kr c n is provided. 
Later, Goldberger and Wise 



23, 



28| introduced a mechanism to provide a potential for 



radion to stabilize the extra dimension. The radion, identitied as the gravitational degree of 
freedom reponsible for the fluctuations of the distance between branes, couples with standard 
model particles with similarity of higgs boson 29M3l| , except that couplings with massless 
gauge bosons could be larger. 

The indications of 125 GeV excess at the LHC has intrigued many interesting discussions 
related with radion 
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371 ] . As shown in [32[, the excess observed at the LHC can be 



explained by a 125 GeV RS radion with a(H)Br(H — > jrf/crBrsM ~ 2.1 and smaller 
values for other channels relative to the corresponding ones in SM. Further phenomological 



discussions are showed in 
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37]. 



The radion couples to SM particles [27J, |28[ as 



where T^ v is energy-momentum tensor for SM particles and A^, = y/6M p ie~ krcn . Thanks to 
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the trace anamoly, this model leads to a larger branching ratio for — >■ gg or 77, relative 
to /ism —> 99 OT 77 in SM. 

In this section, we show that the results of LHC searches for massive graviton have several 
interesting implications for the radion sector in RS model. The first and lightest massive 
Kaluza-Klein(KK) mode of G^ u will couple to SM particles as £j nt = —-^h^T^. As 
is shown above, the couplings of the massive graviton with SM particles are proportional 
to 1/Aq or Xik/Mpi for a fixed Mq. Limits put on Mq for specified k/M p i can then be 
translated to limits on Ac, therefore constraints on due to the relation A^ = \/6Aq. 

Using dijet final states from CMS 2l[ with 1 fb _1 data, we can exclude a RS graviton 
mass below 1 TeV for k/M p i = 0.1 in Fig. [5j A straightforward calculation gives A^, = 
xi ff M - M G = 6.4 TeV. With dilepton final states from ATLAS |22| with 5 fb~\ we show 
in Fig. [6] that a RS graviton mass below 2.2 TeV is excluded at 95% confidence level with 
k/M p i = 0.1, then the corresponding A^ ~ 13.8 TeV. 

A smaller value of A^ then requires a larger k/M p i, although the latter of order 0.1 or less 
is preferred theoretically However, a larger k/M p i means a more stringent constraint on 
Mc because the cross section for the graviton's production at the LHC is proportional to 
(k/M p i) 2 . As shown in Fig. [HJ when k/M p i = 0.3, the limit for Mq is 2.8 TeV, then we have 
A^ = 5.97 TeV. A limit of M G = 3.5 TeV will give A^ = 2.24 TeV for k/M pl ~ 1. Even 
the largest but highly theoretically disfavoured k/M p i ~ 2.88 results in A^ = 0.8 TeV and 
a(H)Br(H -> ^)/aBr SM ~ 1.5. 



V. SUMMARY 

In this work, we discuss the constraints on massive graviton with the latest LHC data. 
In a general framework, we show both dijet and dilepton limits on the two parameters, the 
mass Mq and coupling strengh A G - The limits are applicable to a wide class of models. 
As an illustration, we discuss the implications for RS massive graviton. For k/M p i = 0.1, 
the dilepton search at ATLAS with 5 fb _1 excluded Mq < 2.2 TeV regions. The constraint 
becomes more stringent for larger k/M p i. 

As a byproduct, we show in RS model that constraints on massive graviton can give 
interesting implications on the radion sector. For k/M p i = 0.1, the dijet search at CMS with 
1 fb _1 excluded < 6.4 TeV intervals. And the dilepton search at ATLAS with 5 fb" 1 can 
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exclude A^ < 13.8 TeV regions. For k/M p i ~ 1, the low limit get relaxed to A^, ~ 2.24 TeV. 
These limits have a direct impact on the cross secction a(pp — > <ft) of radion production at 
the LHC since a(pp — > <fi) oc 1/A|. 
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VI. APPENDIX 



A. Decay rate 



This subsection lists the decay rates for graviton to two standard model particles. The 

nn 

complete descriptions and full Feynman rules are refered to [5|, |6[. The decay rate to two 
fermions, G — > f f, is 

where N c is equal to 3 for quarks and 1 for leptons. For vector weak bosons final states, the 
rate is 

r (G -+ WW/ZZ) = 6^j- (1 - Ax v y + + 44) , 

where 5 = 1(|) for W(Z), respectively. For the massless final states, gluon and photon, we 
have 

M 3 

Y{G^gg/ 11 ) = N G 



Here N G = 1(8) for 7(g). Finally, the decay rate to higgs bosons is 
In all the above formulas, Xi = mj/M G , i = f,V, H. 
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B. Bayesian inference 



In this subsection, we briefly outline the statistics used the the context. Based on Bayesian 
inference, given the data or events observed, one can estimate the probability density of the 
parameter s (the signal or cross section, for example) predicted by a theoretical model, 

. . , . £(data|s)7r(s) 

p(s data) = — '/ — . 

Af 

p(s|data) is the posterior probability density function(PDF), £(data|6, s) is the likelihood 
function, tt(s) is the prior PDF, and Af is the normalization constant, 

Af = p(data|s) = J £(data|s)ds. 

Upper limit x can be put on s with 95% confidence level when 

p(s|data)ds = 0.95. 



In a counting experiment, events are the sum of background and singal, the likelihood 
function is actually £(data|6, s) which is given by 

£(data|M = Poisson(#,s) = ^ + g ^ e" (b+s) , 

d\ 

where d is the number of the observed events, b and s are predicted by theory for background 
and signal, respectively. For example, when d = 2, b = 0.92, then the above formalism give 
the upper limit with 95% confidence level s < x = 5.45. Since s is proportional to the cross 
section which is a function of the parameters in the model. Then limit set on s can be 
translated to the parameters. 

For data collected into iV bins, the above formulas can be generalized to 

N AT , 

(Pi 

dA 



£(data|6, s) = JJ Poisson^, Si ) = II ~ t p^ e~ (bi+s °, 



i=l i=l 

where d^b^ s, corresponds the quantities in the 2-th bin. Details and examples can be found 



in [38[ where uncertainty effect is also discussed. 
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